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SUMMARY 

Ravleiah-step helium-purge, annular shaft seals, ^tudi ed for use i n 1 i q- 
uid oxygen turbopumps, generate a used to reduce 

5:re;rp"d\rHe?hanUa, recPPO.opy 

::p'up 

?rIlow airCcrLsed wUh increLed aaial pressure gradient, 

^°hlrthe'’reakrgr "S?edirt?oni^i.aL'a'ma«i™u™'“rrtl,il^ clerrlncptdue 

rSI ' 

for the discrepancy found between measured and predicted flow rates are dis 
cussed. Further testing is recommended. 


INTRODUCTION 

R;i\/1 pi oh-steo helium-purge, annular seals are being studied for use in 

lirliX'e^i^Sa’ncI^-McP resuit in 

reduced leakage. 

Iindpr a contract directed by NASA Lewis Research Center (NAS3-23260) , 

"^d^t1rls5lft"?n^s;rn?nrRarei^"-‘ter:™E^^ 

and fabricated test seals for 20-mm and 50-- shafts and compiled a small 
amount of experimental data for the 50-mm shaft seal. 

Presently an experimental study is underway at NASA Lewis to complete 

crsn! RAsicr^ -.cCicfes^tbe 

ci^C'ificrs\ccr°sccCr"Cs's:yCirci:rsC;crc?iMC iw. 'ihe- other ^code, 

CaefUccaCC‘Cirh%scruCccCe:u:rteHhr:«CurCrfio/Ca;Cs%o^ the 

i“•»;s:»;rsci;cs 

<;pfll that was tested and the analysis that was used in FLOWCAL are p 
ICntCC Cfte pCeCutions. to illustrate the effects various test and 


manufacturing parameters have on the leakage rate, and comparisons of predicted 
and measured flow rates follow. Possible causes for the discrepancies are dis- 
cussed as well as a plan of action to completely validate the use of FLOWCAL 
tor this seal design. 


SYMBOLS 

A annular area, in.- 

C concentric radial clearance, in. 

C 0 discharge coefficient 

D shaft diameter, in. 

e eccentricity, in. 

9c 386 1 bm-i n . /I bf-sec^ 

ft seal dam length, in. 

m mass flow rate, Ibm/sec 

m m/m 

o 

mQ P^C^d + 3/2 e^)A/(24 RT^ jaSl) 

P pressure, psia 

Pa downstream pressure , psia 

^a ^a^^O' riormalized downstream pressure 

Pc Pressure at which choked flow occurs, psia 

^c Pc^Pq- normalized downstream pressure at which choked flow 
occurs 

Pd Ps-Pa. pressure drop across seal, psid 

Pi intermediate pressure, psia 

Pi Pi^Pq’ normalized intermediate pressure 
Pq reference pressure, psia 
Ps inlet stream pressure, psia 
R gas constant, 1 bf- i n . / 1 bm-°R 
Re Reynolds number 
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T temperature, “R 

Ti intermediate temperature, °R 

Ti Ti/T 5 , normalized intermediate temperature 

T 5 upstream or inlet temperature, °R 

V velocity, in. /sec 

X distance, in. 

(5 2RTs (mo/APs)2/gc 

Y specific heat ratio 

e e/c , eccentri ci ty ratio 

y, fluid absolute viscosity, Ibf-sec/in.^ 

p fluid density, Ibm/in.^ 


RAYLEIGH-STEP ANNULAR SEAL 

Pimirp 1 i<; a Dhotoqraph of the 50-mm Rayleigh-step shaft seals and runner 
that «e?rtelted. Each ?ing Is made of Purebon P5N carbon graphite to mimmize 
the mass moment of inertia of the seal, thereby improving the seal s capability 
Jo track sSt perturbations. Four capacitance probes (two diametrically 

Si?ert?y'’''^he Ta^rs'were'y^' apfrt t^k^i U tl?e"e?centrlci t^^e^ 

Arp in-italled A wavy washer separates the two seal rings and holds them 
ao!insr?he housing To ensere the seals' ability to respond to shaft pen r- 
bations the frictional force between the seal rings and the housing is 
mi 7 pd hv a small axial lip (and hence small contact area between the seal g 

an^housing) radially positioned to nearly balance the opposing axial pressure 

forces . 


FLOWCAL 

Theory 

FLOWCAL calculates either isothermal or adiabatic gas flow rates 
an annular clearance with an axial pressure gradient. Eccentricity and inlet 
inertial effects are accounted for. The seal is divided into two regions, 
inlet region where the flow is strictly inertial and a film region where vis 
coisVor?es dominate (see fig. 4). At the inlet, the flow rate is computed 
in a flow nozzle by the following equation; 
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-■ vM " 


^ kj ' - ^ ^ 


( 1 ) 


The critical pressure for choked flow is 

2 


P = 

C 1 Y + 


( 2 ) 


If the normalized pressure downstream of the inlet, Pj, becomes less than the 
critical presure P^. then P. is set equal to ^ ^^d the mass flow rate 

becomes 




1 /2 


2 ^ (y+1 )/2<y-1 ) 


S Vy + 1. 

In the film region, the following equations apply: 


dP 1 dv\ . 

d7 " ^ d^j 


-12 uft 


^ I 


C 


(3) 


(4) 


m = pvA, P = p 

Combining equations (4) and (5) and normalizing yields 


p (y-1)/y 

Iy i /^p <1 +y)''y _ P (1+y)/y\ _ M 

i a 


1 + 


T. 




’ Y 

vpj 


rti 


(5) 


(6) 


The critical pressure at which choking will occur is found by differentiating 

a ' 


with respect to P^, setting di1i/dP^ = 0, and solving for P^. This value of 


P, is the critical pressure, P , where 

d C 


p = 
c 


0*2 T. 

3iti 1 


y/( 1 +yl 


2y p (y-11/y 
i 


( 7 ) 


Again, if P^ < P^, then set P^ = P^. If the flow is isothermal, then 
Y = 1 in the film region equation. 


4 



To compute the flow rate through the seal, a J ^ s 

r;:o'i:rcu?ar"the .als f,o» rates in both the ^ ^ Upon 

*i:'cr?^:rtef u^tii^i: btrt:;e^crse^::rn t^i- ;:t:rn‘:?tS? ^ven 

?Xance r?he maximum number of iterations permitted has been reached. The 
iteration process is an accelerated Newton-Raphson method. 


The following major assumptions are made in this analysis 


(1 ) The annular clearance i s 

(2) The gas is an ideal gas. 

(3) The viscosity and specif! 


small compared to the shaft radius, 
c heat ratio of the gas remain constant. 


Input and Output 

Some of the input required for FLOWCAL is illustrated in figure J^e 
fluid conditions and properties needed are the upstream pressure and tempe 

r>f thp discharae coefficient is between 0.6 and 1.0. The program has the 
option to calculate adiabatic or isothermal flow. The program output ^ 

?pmnp?ature pressure and Reynolds number; choked flow pressure, flow rate, and 
fS er^r* The^now er?or is the difference between the flow rates calculated 
from the inlet region equation and the film region equation. 


ANALYTICAL RESULTS AND DISCUSSION 

Parametric predictions were made to illustrate how the flow rate is 
affected by operating conditions and to bound the predicted flow rate range by 
michininq tolfran es and reasonable levels of measurement uncertainty. Predic- 
Hon were mad^for both the 50-mm and 20-mm diam shaft sizes. However only 

fied on the figures. 

The baseline cases for adiabatic and isothermal flows are shown in fig- 
ure 5 The volumetric flow rate (leakage rate) in standard cubicfeet per 

minute (SCFH). which it directly proportlonai to mass ^ “o2 

fiinrtinn of the oressure drop across the seal P(j. As expectea, tne 
. i^c?easL with a^^ncrease in Pd- Leakage rates for adiaba ic f ow are 
greater than for isothermal flow and the difference increases wi me 

Pd- 

The effect of the downstream pressure on flow >^®te is shown in figure 6. 
PrpHirtions were made for three downstream pressures: 14.7, 4b. u, ana 

75 0 psia. Note that the upstream pressure is the sum of the 
sure and the pressure drop across the seal. Flow rates are greater for t 


a 

rate 
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higher downstream pressures because the fluid density is greater. At a pres- 
sure drop across the seal of 100 psid or more, the slopes of all three curves 
are nearly equal . 

Figure 7 shows the leakage rate as a function of inlet stream temperature 
for P(j equal to 40.0 psid. The correct viscosity for each inlet temperature 
was used. Theflow rate decreases as the inlet temperature increases because 
the fluid density is inversely proportional to temperature as stated in the 
ideal gas law. Superimposed on the linear effect of density change is the 
nonlinear effect of viscosity, which increases with temperature and further 
reduces the leakage rate. 

Figure 8 is a plot of predicted flow rate with respect to eccentricity 
ratio e. The leakage rate is only slightly affected by an increase in eccen- 
tricity ratio as long as c £ 0.2. For e > 0.2 the flow rate is dramatically 
increased. This result is consistent with the experimentally verified analysis 
done by Tao and Donovan (ref. 1). The ratio of the fully eccentric flow rate 
to the concentric flow rate for the pressure drop across the seal of 5.0 and 
115.0 psid are 2.466 and 1.564, respectively. Calculation of the Reynolds 
number shows that at P(j = 5.0 psid the flow is laminar (Re = 33.5) and that at 
P(j = 115.0 psid the flow is turbulent (Re = 4415). In reference 1, analysis 
and experimental data show that the fully eccentric to concentric flow rate 
ratio should be approximately 2.5 for laminar flow and between 1.2 and 1.4 for 
turbulent flow. This supports the predictions made with FLOWCAL. 

In an effort to bound the range of flow rates to be measured, predictions 
were made for the extreme cases of centrifugal and thermal growths, machining 
tolerances, and a ±2 percent uncertainty in the clearance measurements. All of 
these items change the range of clearances at which the seal operates. Since 
the flow rate is a function of the clearance cubed, minor changes in the clear- 
ance can have significant results. 

The amount of centrifugal growth of the 50-mm diam runner for the 
Rayleigh-step, helium-purge seal may be found in reference 2. Predicted flow 
rates for five different shaft speeds are shown in figure 9. As expected, flow 
rates are lower for the higher shaft speeds because the seal clearance has been 
reduced by centrifugal growth of the shaft. 

In order to band the range of flow rates due to thermal growth effects of 
the seal and runner, the two extreme cases were considered. The maximum clear- 
ance due to the thermal effects is obtained when the shaft is at the tempera- 
ture of liquid nitrogen (-320 °F) and the seal ring is at the temperature of 
ambient helium (70 °F). (Liquid nitrogen hydrostatic bearings are used in the 
seal tester to support the shaft; hence the shaft gets extremely cold.) The 
minimum clearance occurs when the shaft and seal ring both operate at ambient 
(70 °F) conditions. The flow rates predicted for these two conditions are 
shown in figure 10. 

Another way to bound the flow rates is to make predictions based on the 
machine tolerances allowed in the fabrication of the shaft runner and seal 
ring. This is illustrated in figure 11. At a 200.0 psid pressure drop across 
the seal, the flowrate could be between 27.5 and 15.9 SCFM, which represents a 
difference of 42.2 percent. 
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In determining how the flow rate is affected by error in the clearance 
measurement, predictions were made for the baseline clearance ±2 perce 
These predictions are shown in figure 12. Compared to thermal effects or 
machining tolerances, the range of flow rates is bounded much more close y, 
with a maximum variance of only 1.75 SCFM. 


COMPARISON BETWEEN PREDICTED AND MEASURED FLOW RATES 

Flow rates were predicted for the 50-mm Rayl ei gh-step , helium-purge seal 
by using clearance, eccentricity, and pressure measurements made during test- 
inn a-; innut to the analysis. The inlet temperature was assumed to be 70 F. 
Tese predicted flow ratL are plotted against the measured flow |;jtes in fig- 
ure 13 If predictions and measurements were in complete agreement. 

<;hould oroduce a straight line with a slope of 1.0 that passes through the on 
should produce a sc a g measured flow rates are approxi- 

mately three to’seven times greater than predicted, and the data follows the 
rurve of ^second-degree polynomial. This large discrepancy forces one to 
Question the validity of the analytical model for this application^ To verify 
that ^e analys does accurately predict the gas flow rate through an annu ar 
Teal !ith an Lial pressure gradient, data independent of this experiment was 

sought . 

Fiqure 14 shows data taken by Dike (ref. 3) with analytical 
<;nnprimLsed. A discharge coefficient of 0.6 was used. Agreement between the 
analysis and this data is very good. Also, data published by re . 

chni . mil^r aarLment. (See fig. 15.) In this case, discharge coefficients 
of 1 0 and 0 8 were used in predicting flow rates. To determine the reasona 
bleness of using the selected discharge coefficients, predictions of the coef 
ficients were made by using theoretical and empirical equations for annular 
or f i e Tref 5) For the configuration of Dike's work (ref^ 3), the dis- 
chlrqrcoefficient is predicted to be between 0.68 and 0.55 depending on if 
tSrorf ?fi ha p-edged or square-edge and thick, respectively. Thedis- 
hL^e oeffm of choice, 0.6. falls in this range. For the configure ion 
used in Nelson's (ref. 4) work, the discharge coefficient is predicted to be 
between 0 67 and 0.71, which is lower than what is needed to match the expe 
mental data in figure 15. This discrepancy may be due to the ^°^ 9 er length of 

r eaf and the om. "itTlJerSi^rt" ^^^t ThrflS^s^^ 

Sed „as not choKed in the 50-™ 

Rayleigh-step, helium-purge, seal tests. 

In considering the possible explanations for the discrepancy between the 
predicted and measured flow rates for the 50-mm Rayl ei gh-step . helium purge 
seal several possibilities come to mind. 

First the unchoked flow portion of the analysis could be ignoring phenom- 
ena that significantly affect the flow rate, although the assumptions, equa- 
tions and trends found in parametric studies are reasonable. Unfortunately, 
no experimental data for unchoked flow through an annular clearance found 
so that the analysis could be completely verified at this time. It is appa 
ent however, that the analysis does accurately predict flow rates for choked 
flow This may be because the flow nozzle equation at the inlet dominates the 
film region equation when the flow is choked, which makes the flow rate less 
sensitive to the clearance measurement. 
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Second, there could be errors in the experimental data. Error in the 
measured clearance is one parameter that would have significant effects because 
flow rate is a function of the clearance cubed. The smallest possible clear- 
ance would occur if the shaft and seal ring were at ambient temperature and at 
maximum speed. The largest clearance would occur if the shaft temperature was 
that of liquid nitrogen and the seal ring remained at ambient temperature. 
Although it is not probable that these exact conditions did occur, they do pro- 
vide a range within which the measured clearances should be. All the clearance 
measurements were between these two extremes of 0.00018 to 0.002934 in. 

A third possible explanation for the discrepancy between measured and pre- 
dicted flow rates may be that a second flow path existed between the housing 
and the seal ring. This could be due to imperfections in either the seal or 
the housing surface or both. Also, the carbon seal ring could be exhibiting 
an unstable geometry; twisting when exposed to cryogenic temperatures and 
large thermal gradients. In addition, MTI observed an inboard seal ring that 
developed an elliptical orbit in phase with the shaft orbit. They speculated 
that the frictional resistance of the seal ring was minimized by the existence 
of a fluid film between the ring and the end wall (ref. 2). 


CONCLUDING REMARKS 

Having no data to verify the unchoked flow analysis, and knowing there is 
the possibility of a secondary flow path and that small errors in the clearance 
measurement have great effect on the predicted unchoked flow rate, I recommend 
further testing be done. In addition to closely checking clearance and flow 
rate measurements, this testing should entail inspection of the axial lip and 
housing sealing surfaces for imperfections. Testing should be done at both 
room and cryogenic temperatures while particularly looking for twisting of the 
carbon seal at low temperatures. Also, a static test of the 50-mm, helium- 
purge seal and the seal runner configured in an apparatus that insures no sec- 
ondary leakage at the axial lip should be conducted. Flow rates and clearances 
for several axial pressure gradients should be measured. If flow rates of this 
configuration were significantly less than those previously measured, I would 
conclude a secondary flow path did exist in previous testing. Comparisons of 
predicted and static test flow rates should then be made to verify the code. 

It is important to measure flow rates for both choked and unchoked flow to com- 
pletely verify the analysis. By following these steps, the cause of the dis- 
crepancy between the analysis and measured data should be revealed. 


SUMMARY OF RESULTS 

Analytical data showed that the gas flow rate through an annular clearance 
with an axial pressure gradient increased with increased pressure drop across 
the seal. Flow rate also increased with higher downstream pressure because the 
density is higher; flowrate is greater for adiabatic flow than for isothermal 
flow. Flow rates decreased with elevated inlet temperature because of lower 
density and higher viscosity. Also, predictions found flow rates increased 
with eccentricity ratio, which is consistent with previous research. Study 
of those things that affect the radial clearance, such as centrifugal growth, 
thermal growth, machining tolerances, and clearance measurement uncertainty, 
placed wide boundaries on the expected flow rates. The widest boundaries were 
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due to thermal growth effects^ Predicted flow rates underestimated measured 
flow rates for a 50-mm Rayleigh-step, helium-purge seal by three to seven 
tSes However the analysis accurately predicted the flow rates measured and 
JZ?ted TotlJe. e? al ^ and Nelson, et al . for choked gas flow through annular 
seals. Further testing is recommended. 


TABLE I. - BASELINE INPUT FOR PREDICTED FLOWRATES 
THROUGH THE 50-mm RAYLEIGH-STEP. ANNULAR SHAFT SEAL 


[Isothermal flow] 


Downstream pressure, psia 
Upstream temperature, “F . 

Seal length, in 

Shaft diameter, in. . . 
Radial clearance, in. 
Eccentricity ratio . . . 
Discharge coefficient 
Viscosity, Reyns .... 
Ratio of spec i f i c heats 
Molecular weight .... 


. . . I ~r . > 

. . . 70.0 
. . 0.077 

. . 1.9684 
. 0.000669 
... 0.0 
... 1.0 
2.8855x10-9 
. . . 1.66 
. . 4.003 
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C-85-4966 


Figure 1 - Rayleigh-step helium-purge annular seal pair with 50-mm runner. Ring is made of Purebon 
P5N carbon graphite and runner is Inconel 718. 


/ — CIRCUMFERENTIAL GROOVE 


SHAFT 



ORIGINAL PAGE IS 
OF POOR QUALITY 


Figure 2. - Inner surface of the 50-mm Rayleigh-slep, annular seal ring. (AN 
linear dimensions are in Inches.) Inside and outside seal ring diameters are 
1.9694, and 2.66 inches, resp>eclively. 





figure 5. - Effect of isothermal and adiabatic flow on predicted 
leakage rate as a function of pressure frop across the 50- 
mm nayleigh-step. helium purge, annular seal. 



Figure 6. - Effect of downstream pressure on predicted leakage 
rate as a function of pressure drop across the 50- mm 
Rayieigh-step, helium purge, annular seal. 
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INLET TEMPERATURE. 

Figure 7. - Effect of inlet stream temperature on predicted 
leakage rate of 50-mm Rayleigh-step. helium purge, annular 
seal for pressure drop of 40.0 psld across seal 



Figure 8. - Effect of eccentricity ratio and predicted leakage 
rate of 50-mm Rayleigh-step, helium-purge, annular seal for 
pressure drops across seal of 5.0 and 1 15.0 psid. 
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SPEED. 

rpm 

SHAFT 

DIAMETER. 

In. 

RADIAL 
CL LA RANG! 

in 

o 

0 

1.9684 

6.69 X 1o"^ 

□ 

30000 

1-96856 

5.89 X lO"'^ 

o 

45000 

1 .968776 

4.01 X 1o"^ 

A 

60000 

1.96906 

3.39 X 10 ^ 

V 

70000 

1 .9693 

2.19 X 10 



20 40 60 00 100 1 20 140 160 180 

PRESSURE DROP ACROSS SEAL, psid 


Figure 9. - Effect of centrifugal growth caused by shaft speed 
on predicted leakage rate as a function of pressure drop 
across the 50-mm Rayleigh-step, helium purge, annular seal, 


DIAMETE R. CLEARANCE . 



40 60 80 100 120 140 160 

PRESSURE DROP ACROSS SEAL, psid 
Figure 10. - Maximum and minimum predicted leakage rates 
caused by thermal growth changing the clearance of the 50- 
mm Rayleigh-step. helium purge, annular seal, 
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LEAKAGE RATE, SCFM PREDICTED LEAKAGE RATE, SFCM 
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DIAMETER - 


0^ 60^ 80 100 120 140 160 180 200 

PRESSURE DROP ACROSS SEAL, psid 
FiQure 11- Maximum and minimum predicted leakage rates 
caused by machining tolerances tor 50 mm Rayleigh step, 
helium purge, annular seal 
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PRESSURE DROP ACROSS SEAL, psid 

Fiaure 12 - E fleet of a ±2 percent uncertainty in clearance meas- 
urement on predicted leakage rate for 50 mm Rayletgh-step, 
helium-purge, annular seal. 
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0 2 4 6 8 10 12 14 16 18 2C 

MEASURED LEAKAGE HATE, SCFM 
Fiaure 13 - Predicted leakage rates as a function of measured 
leakage rates for 50-mm Rayleigh step, helium-purge, 
annular seal, compared to perfect prediction line. 


) 400 800 1200 1600 200C 

SUPPLY PRESSURE, psia 

Fiaure 14 - Comparison of predicted leakage rate with leak- 
age rate measured by Oike et al (ret. 3) for a Rayleigh- step, 
gaseous-hydrogen, annular seal. Inlet ~ ^^^,.8' 

downstream pressure = 14.7 psia, viscosity - 042 
Reyns, specific heat ratio = 1 .39. seal length = 0.28 in., sha t 
diam - 1 97 in., radial clearance = 2,26 x 10 in., eccentricity 
ratio = 0.0. discharge coefticient = 0.6. molecular weight = 

2.0, adiabatic flow. 


ANALYSIS Cf3 - 

NO PRE ROTATION 
POSITIVE PRE ROTATION 
NEGATIVE PREROTATION 


I I .... I I ‘ ‘ 

, 20 40 60 80 101 
PRESSURE DROP ACROSS SEAL, psid 
Figure 15 - Comparison of predicted and measured leakage 
rates for sea! tested at Texas A&M University (ref. 4). Inlet 
temperature = 70 '^F. downstream pressure = 14.7 psia, 
viscosity = 2.755 x 10'9 Reyns, specific heat ratio ^1 4. 
molecular weight - 28.9, sea! length = 2.0 in,, shaft diam^ 

6 0 in., radial clearance = 0.029 in,, eccentricity ratio - 0 0, 
isothermal flow. Positive prerotation means the inlet flow 
had a circumferential velocity in the direction of shaft 
rotation. Negative prerotation is against the direction of shatt 
rotation. 


APPENDIX - FLOW RATE PREDICTIONS FOR 20-mm RAYLEIGH-STEP 
HELIUM-PURGE, ANNULAR SEAL 

TABLE II. - BASELINE INPUT FOR PREDICTION OF FLOW RATES 
THROUGH THE 20-mm DIAMETER, RAYLEIGH-STEP, ANNULAR 

SHAFT SEAL 

[Isothermal flow] 


Downstream pressure, psia 14. 7 

Upstream temperature, °F 70. o 

Seal length, in 0.077 

Shaft diameter, in 0.78745 

Radial clearance, in 0.000325 

Eccentricity ratio 0.0 

Discharge coefficient 1.0 

Viscosity, Reyns 2.8855x10-9 

Ratio of specific heats 1.66 

Molecular weight 4.003 



Tigurt' 1fi. rffoci of isothermal and adiabatic flow on predicted 
leakage rate as a function of pressure drop across the 20-mm 
Mayleigh step, helium-purge, annular seal. 




PREDICTED LEAKAGE RATE, SCFM PREDICTED LEAKAGE RATE, SCFM 



Figure 17. ^ Effect of downstream pressure on predicted 
leakage rate as a function of pressure drop across 20 mm 
Rayleigh-step. helium-purge, annular seal. 



INLET TEMPERATURE, '’F 

Flqure 18 - Effect of inlet stream temperature on predicted 
leakage rate of 20-mm Rayleigh-step. helium-purge annular 
seal with a pressure drop across the seal of 40.0 psid 




FCCLNTRICITY RAIIO 


Figure 19. Effect of eccentricity ratio on predicted leakage 
rate of 20-mm Rayleigh-step. helium purge, annular sea! tor 
pressure drops across seal of 5.0 and 85.0 psid. 
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Figure 20. - Maximum and minimum predicted leakage rates 
caused by thermal growth changing the clearance in the 20 
mm Rayleigh-step, helium purge, annular seal. 
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Figure 21. - Maximum and minimum predicted leakage rales 
due to machining tolerances for 20 mm Rayleigh step, 
helium-purge, annular seal. 
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